Theory of the Lightly Doped Mott insulator 
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A theory for the Hubbard model appropriate in the limit of large U /t, small doping away from 
half-filling and short-ranged antiferromagnetic spin correlations is presented. Despite the absence of 
any broken symmetry the Fermi surface takes the form of elliptical hole pockets centered near (-f, §) 
with a volume proportional to the hole concentration. Short range antiferromagnetic correlations 
render the nearest neighbor hopping almost ineffective so that only second or third nearest neighbor 
hopping contributes appreciably to the dispersion relation. 
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I. INTRODUCTION 

The existence and shape of the Fermi surface and 
its change with doping may be one of the central is- 
sues in the physics of cuprate superconductors. Quan- 
tum oscillations consistent with Fermi liquid behaviour 
and a pocket-like Fermi surface have been observed 
in the underdoped compounds YBa2Cu306.5[i|-13| and 
YBa2Cu408[5|, Q. Overdoped Tl2Ba2Cu06-i-5 on the 
other hand shows quantum oscillations as well but a 
'large' Fermi surfaceJT] consistent with band structure 
calculations. This indicates a change of the Fermi surface 
volume around optimal doping. Similarly, the electron 
doped compound Nd2_a:Ce3:Cu04 shows a pocket-like 
Fermi surface for electron concentrations up to (5 = 0.16 
and which changes abruptly to a large Fermi surface for 
5 — 0.17[8J. The transition in Nd2-a:Cej;Cu04 thus oc- 
curs in an overdoped compound and therefore is unre- 
lated to antiferromagnetic ordering. 
The validity of a Fermi liquid description in underdoped 
cuprates is incompatible with the 'Fermi arc' picture 
which is frequently invoked to describe the absence of 
a large Fermi surface in angle resolved photoemission 
spectroscopy (ARPES)[9[. The Fermi surface of a Fermi 
liquid is a constant energy contour of the quasiparticle 
dispersion and thus necessarily a closed curve in k-space. 
ARPES experiments on insulating antiferromagnets like 
Sr2Cu202Cl2[ll and Ca2Cu02Cl2[ll| have shown that 
the valence band is consistent with next-nearest hopping 
- as in an antiferromagnet - with maximum close to (^, ^) 
and that the part of the quasiparticle band facing (tt, tt) 
has very small spectral weight. Assuming that the ef- 
fect of doping mainly consists in the chemical potential 
cutting into this quasiparticle band the Fermi surface 
would take the form of elliptical hole pockets and the 
'Fermi arcs' would simply be the portions of the pocket 
with large spectral weight. This has in fact been con- 
firmed by the recent ARPES experiments on underdoped 
Bi2(Sr2--2,La2,)Cu06 by Meng et aZ.fT^]. 
As for the quantum oscillation experiments the electron- 
like nature of the carriers suggested by the sign of both 
Hall constant J13J and thermopower'14J is incompatible 



with a straightforward interpretation in terms of hole 
pockets. Rather, the strong temperature dependence 
of the Hall constant [l5| and the enhancement of low 
frequency magnetic excitations in a magnetic field [l6| 
suggest a complicated and as yet not understood re- 
construction process to take place in YBa2Cu306.5 and 
YBa2Cu408. Concerning the thermopower it has also 
been pointed out that in the cuprates there may be no 
obvious correspondence between the sign of the ther- 
mopower and the Fermi surface geometry 17 . 
In the present manuscript we investigate the point of view 
that hole pockets are a generic property of a lightly doped 
Mott insulator where the bulk of electrons continues to be 
localized as in the insulator and the mobile carriers corre- 
spond to the doped holes. The localized electrons retain 
only their spin degrees of freedom and do not contribute 
to the volume of the Fermi surface which leads to a Fermi 
surface with a volume proportional to the hole concen- 
tration irrespective of any broken symmetry. In fact, no 
experimental evidence for any staggered order parame- 
ter which would explain hole pockets by backfolding of a 
large Fermi surface has been found so far. Moreover this 
picture - a single mobile hole interacting with spin exci- 
tations - is the underlying one for all successful theories 
of the ARPES spectra of insulating compounds |18l - [26| . 
Further motivation for the present work comes from ex- 
act diagonalization studies of the t-J model. These 
show that the Fermi surface at hole dopings < 15% 
takes the form of hole pockets |27h29| . that the quasipar- 
ticles have the character of strongly renornialized spin 
polarons throughout this doping range[30-32] and that 
the low energy spectrum at these doping levels can be 
described as a Fermi liquid of spin 1/2 quasiparticles 
corresponding to the doped holes [33]. A comparison of 
the dynamical spin and density correlation function at 
low {5 < 15%) and intermediate ((5 = 30 - 50%) hole 
doping moreover indicates that around optimal doping a 
phase transition takes place. In the underdoped regime 
spin and density correlation function differ strongly, 
with magnon-like spin excitations and extended incoher- 
ent continua in the density correlation function 34 |35| 
which can be explained quantitatively by a calculation in 



the spin-polaron formalism [36j. At higher doping, spin 
and density correlation function become more and more 
similar and both approach the self-convolution of the 
single-particle Green's function, whereby deviations from 
the self-convolution form can be explained as particle- 
hole excitations across a free electron- like ('large') Fermi 
surface [37[. This rough picture would be similar to the 
present experimental situation for cuprate superconduc- 
tors. 

Here we present a theory for the underdoped phase. We 
study the 2 dimensional (2D) Hubbard model 

H = Ht+Hu 






(1) 



where the nearest neighbor hopping tio = 1 and we fix 
U/tio = 8. The results depend only weakly on U/t 
as long as this is sufficiently large. In addition we use 
til = —0.1 tio and ^21 = —tii/2. These values of in 
and ^20 are smaller than the generally accepted ones for 
cuprate superconductors - this will be discussed below. 
In setting up a theory we have the following picture in 
mind: at half-filling - i.e. the Mott insulator - the elec- 
trons are localized and retain only their spin degrees 
of freedom. The hopping term creates charge fluctua- 
tions i.e. holes and double occupancies which we con- 
sider as spin-i Fermions. The excitation spectrum of 
these Fermions has the well-known Hubbard gap of or- 
der U which exsists irrespective of any kind of order or 
broken symmetry. In section II we set up the Hamilto- 
nian for these charge fluctuations. Since we really want 
to study the doped system where apparently no broken 
symmetry exists we thereby consider a hypothetical insu- 
lating phase with no long range antiferromagnetic order 
but short ranged antiferromagnetic spin correlations i.e. 
the 'spin liquid'. The 2D Hubbard model has no bro- 
ken symmetry at finite temperature so we believe it is 
quite reasonable to study the charge fluctuations in such 
a disordered phase. We then assume that for low doping 
S this picture remains applicable, that means the holes 
created by doping have the same character as the holes 
created by charge fluctuations at half-fllling. Whether 
this assumption is justifled is a question to be answered 
by experiment, but we believe that the recent experi- 
mental results lend some support to this idea. Section 
III gives a summary and conclusions. 



II. EFFECTIVE HAMILTONIAN FOR CHARGE 
FLUCTUATIONS 

The basic idea of the calculation is the (approximate) 
diagonalization of the Hubbard Hamiltonian in a suitably 
chosen truncated Hilbert space. As a starting point for 
constructing the truncated Hilbert space we consider a 



state 1^0) which has exactly one electron/site, is invari- 
ant under point group operations, has momentum zero 
and is a spin singlet. These are the quantum numbers 
of a vacuum state and indeed |^o) will play the role of a 
vacuum state, i.e. a state containing neither charge nor 
spin fluctuations. The only property of this state which 
will enter the calculation is the static spin correlation 
function 



Xij 



(vI/o|S,.S,|*o) 



(2) 



We consider Xij as given and do not attempt to compute 
it. We assume it to be antiferromagnetic and of short 
range i.e. 



where Q 
straint 



X,^. = Co e'Q-(^-^.)e-^T^ (3) 

(7r,7r). Moreover Xij has to obey the con- 






(4) 



which follows from \'^o) being a singlet. In practice we 
assume that ([3]) holds only for more distant than nearest 
neighbors and take the nearest neighbor spin correlation 
Xio > —0.33 in 2D 38] as a flrst free parameter. Next 
we choose ( and adjust Co for the longer range part of 
X so as to fulflll (U). The results for the quasiparticle 
dispersion turn out to be almost independent of C and 
only weakly dependent on xio- 

Having specifled the 'spin background' j^^o) we introduce 
the basis states of the truncated Hilbert space. Using 
the familiar Hubbard operators dj ^ 

c]^(l — nj g) they take the form 



cL'^i 



and c 



N^ 



Nu 



2w.+A^.)/2 j|dX^^^ n 



l*n 



(5) 



These states have double occupancies and holes at spec- 
ified positions and we treat these holes and double occu- 
pancies as weakly interacting spin-^ Fermions, which is 
the key approximation of the theory. Fermions are the 
only meanigful description for these particles because the 
Hubbard operators at different sites anticommute. Since 

i'i'old.Jlj^o) = (*o|cl,.c,,J*o) - i the states m 
are approximately normalized if the average distance be- 
tween the holes and double occupancies is larger than the 
spin correlation length ( (see below). This condition is 
satisfied in the limit of large U, small doping and short 
spin correlation length, which is the case of interest. 
A key problem in setting up a theory for the charge fluc- 
tuations in the Hubbard model is the peculiar nature of 
holes and double occupancies. Both a hole and double 
occupany are spinless objects. Despite this, for example 
the states d|-.|\E'o) and (T-^l'^o) are orthogonal. In fact, 
acting with d- ^ implies a projection onto the component 



of l^o) which has a a electron on site i. While the newly 
created double occupancy is a spinless object, the infor- 
mation about the spin of the added electron therefore is 
'stored' in the 'spin background' and it is in a sense stored 
within a spatial region of extend C around the respective 
double occupany/hole. This can be seen by considering 
expressions like 






1 



;Xii 



1 



1 

4 + 3 
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Xij 



*ol4;'^j-.i^],T^^tl*") 



->Xij 



(6) 



and similar ones. If the spin correlation function x were 
zero in these and similar expressions the states ^ would 
indeed be normalized. In the following we neglect the 
corrections due to a finite x in overlap matrix elements 
such as ([5]). This is probably the most problematic ap- 
proximation in the present theory and induces some in- 
accuracies - as will be discussed below. Once we neglect 
the overlap at short distances, however, the states ([5]) are 
indeed normalized. 

Finally, the holes and double occupancies have to obey a 
hard-core constraint i.e. there may be at most one parti- 
cle per site. In the case of small (5/large C//i, however, the 
density of these particles is small so that it is probably a 
good approximation to neglect the hard core constraint, 
see Appendix A for a comparison of the present theory 
with linear spin wave theory where the hard core con- 
straint between magnons is neglected as well. 
The procedure to be applied then is quite simple: the 
states ([5]) are represented by states of fictitious Fcrmionic 
spin-i quasiparticles 



N„ 



N^ 



n4..„ n^i^jo)- 



(7) 



i/=i 



M=l 



This means we have hole-like quasiparticles h\ ^ and 

the double occupancy-like quasiparticles d\ ^ . All op- 
erators in the quasiparticle Hilbert space are defined by 
demanding that they their matrix elements between the 
states ([7]) are identical to those of the physical operators 
between the corresponding states ([5]). 
We now use this procedure to set up the quasiparticle 
Hamiltonian. One has 



1 



(*0|4,d,,, i?t |*0) = -UJ{--X^l) 
(vI/o|4. fft 4^l*0> = -UA\+X^,) 
(*0|C,V ilt €,,,1*0) = UA\+X^j) (8) 

These matrix elements describe the pair creation of a 
hole/double occupancy, the propagation of a double oc- 



cupancy and the propagation of a hole. Thereby modifi- 
cations due to nearby additional particles are neglected 
but again this will be reasonable for small particle density 
and short spin correlation length. Denoting 



ij ~ ''ijin ^Xij) 



y = UAl^+^Xr,) (9) 



the Hamiltonian governing the quasiparticles therefore is 

hL + h,^d. 



E^.(4 



(T'"jcr ' '"JO- 



.0 



«J,o" 



+[/^4,d,,. 



(10) 



The last term takes into account the fact that each double 

occupancy increases the energy by L'^.The extra factor of 

2 in ([S]) as compared to ([5]) takes into account the factor 

of2W^+^-)/Mn(l5l). 

The Hamiltonian ([TU]) is solved by the transformation 



li,^,a = -«k 4.,T + "k ^-kff 

and we obtain the energies 

£;±(k) = e~k + ^±Wk 



(11) 



W^k = \l{^? + vi 



(12) 



and the coefficients 



Mk 



'Wk+i 
2W^k 



'M^k 



Wk 



2W^k 



sign (Vic) 



(13) 



Since Ck,cr = rfk,CT + Ck,cr the representation of the electron 
annihilation operator in the quasiparticle Hilbert space 
becomes 



Ck,a -^ ^(4,^ + /l-k,ff) 



(14) 



where the factor of l/-\/2 again is due to the prefactor in 
([5]). The spectral weight of the two bands therefore is 



Z±(k) = -(uk±«k)' 
2 V VKk 



(15) 



If we set Xij = both e^ and — Vk reduce to ek/2, with 
eic the noninteracting band energy and we obtain 



£;±(k) 



z±(k) 




(16) 



which is identical to the result of the Hubbard-I 
approximation [39J at half-filling. The present theory thus 
may be viewed as an extension of the Hubbard-I approxi- 
mation to take into account the effect of finite spin corre- 
lations. The spin correlations, however, do have a drastic 
effect: for xw = —0.25 the nearest neighbor hopping ma- 
trix element tio is zero. Here one has to bear in mind that 
the value of xio in the ground state of the Heisenberg an- 
tiferromagnet on the 2D square lattice with nn exchange 
only is « — 0.33[38j - a slight reduction of the spin cor- 
relations due to hole doping and longer range exchange 
may well produce a value very close to —0.25. For small 
tio, however, the quasiparticle dispersion is dominated 
by next-nearest neighbor hopping and the quasiparticle 
dispersion is 'almost antiferromagnetic' and thus quite 
different from the Hubbard-I approximation. 
The second major difference between the present the- 
ory and the Hubbard-I and similar approximations is the 
way in which electrons are counted. In the quasiparticle 
Hilbert space the operator of electron number obviously 



N^=N- 



E(< 



^,ad^,a - KaKa) 



(17) 



because the 'spin background' |^o) contributes N elec- 
trons (with TV the number of sites), each double occu- 
pancy increases the number of electrons by one and each 
hole decreases the number of electrons by one. When 
applied to a quasiparticle state of the type ([7]) the op- 
erator p?|) therefore gives the same electron number as 
the physical electron operator applied to the correspond- 
ing Hubbard model state (O and this is the prescription 
how operators in the quasiparticle Hilbert space are to be 
constructed. After transformation to the 7's and treating 
these as noninteracting Fermions we obtain 



N, 



Ei 



,t 



7i^,+,.7k,+.. + 7iL,-,a7k,-,. ) - N. (18) 



At half- filling, N^ = N, the lower of the two bands is 
completely filled, the upper one completely empty, which 
agrees with the Hubbard-I approximation. As the sys- 
tems is doped away from half- filling, however, the Fermi 
surface volume has a volume which is strictly propor- 
tional to the number of doped holes - i.e. one has hole 
pockets with a total volume of 5/2. These do not occur 
as a consequence of backfolding the Brillouin zone due to 
any kind of broken symmetry. 

Pockets with a volume of S/2 are different from the 
Hubbard-I approximation and the reason is that there 



one uses a different way to count electrons, namely the 
integrated photoemission weight. Treating the d and h 
as ordinary Fermion operators, and using ^T^ we obtain 
the integrated spectral weight as 



k,CT 



4E(4> 



-kff ' "-kCT^lcCT- 



(19) 



which differs from (fT7|) . There are several reasons for this 
discrepancy: first of all we have 



E< 



4,>ik^ 



= E(4>i 



(20) 



The latter expectation value, however, should be zero, 
because a hole and a double occupancy cannot occupy 
the same site. Even with this term omitted, however, 
there is an extra factor of 1/2 and reason is a more fun- 
damental one, namely the restriction of the Hilbert space. 
In fact it is easy to see that the spectral weight sum-rule 
cannot be applied to an approximation like the present 
one where the spectral weight is artificially concentrated 
in a single band. Let us consider a Fermi liquid with 
quasiparticle weight Z < I and assume that one electron 
is removed at the Fermi energy. This implies that one 
momentum/spin crosses from the occupied to the unoc- 
cupied side of the Fermi energy which decreases the inte- 
grated photoemission weight by Z. The remaining weight 
of 1 — Z therefore must disappear from the high-energy 
part of the photoemission spectrum and reappear in the 
high-energy part of the inverse photoemission spectrum. 
This is in fact exactly what is seen in exact diagonal- 
ization studies of the t-J model: there on has a quasi- 
particle band of width « 2J and quasiparticle weight 
Z « 0.1 — 0.5 while most of the spectral weight resides 
in extended incoherent continuaj4Cl |4l|. Upon doping 
the quasiparticle peak at the top of the band crosses the 
chemical potential while simultaneously spectral weight 
is removed from the incoherent part of the spectrum at 
energies of order t below the Fermi energy. This weight 
reappears - in the form of multi-magnon excitations - at 
energies of order J above the chemical potential and near 
(tt, tt) [30| (another way is realized in SDW mean- field the- 
ory for the Hubbard model, where the weight Z crosses 
at k and the remaining weight 1 — Zatk-I-Q). In a 
theory like the present one where all spectral weight is 
concentrated in one quasiparticle band, however, these 
high-energy parts do not exist and it is therefore impos- 
sible to maintain the spectral weight sum rule. In fact, 
for each electron removed from the system, Z~^ momenta 
would have to cross from photoemission to inverse pho- 
toemission to maintain the sum-rule and it therefore is 
misleading to use the integrated photoemisssion weight 
to determine the Fermi surface volume. We have to ac- 
cept that the spectral weight sum rule will not be fulfilled 
as a consequence of the restriction of the Hilbert space 



and the lack of the incoherent part of the spectra (this 
may be remedied at least partially by including the cou- 
pling to spin excitations so as to reproduce the incoherent 
continua). Instead, the correct expression for the particle 
number is given by (|17l) . 

In many calculations where the Hubbard-I approximation 
- or any other approximation involving Hubbard opera- 
tors - is applied to the doped case, the analogue of P^ 
is used to calculate the electron number. This leads to 
the peculiar 'fractional' dependence of the Fermi surface 
volume on hole doping because of the fractional num- 
ber Z^^ of momenta needed to fulfill the spectral weight 
sum-rule. 

We proceed with a qualitative discussion of the band dis- 
persion. Since we expect the theory to be valid only for 
t ^ U we thereby expand in terms of t/U and keep only 
the lowest nonvanishing order. Moreover we note that 
there are two additional small parameters. For the near- 
est neighbor spin correlation function xio close to —4, 
the effective nearest neighbor hopping tio is small. The 
value of xio in the GS of the 2D Heisenberg antiferro- 
magnet with nearest-neighbor exchange is — 0.33[38], so 
that a slight reduction due to doping/longer range hop- 
ping may well produce a value which is very close to —7. 
Moreover, in the Cu02 plane the (1,1) and (2,0) hop- 
ping integrals fulfill ^11/^20 ~ ^2. This equation holds 
if the underlying mechanism for these matrix elements 
is hopping of a Zhang-Rice singlet via Cu 4s orbitals. If 
the spin correlations xn and X20 do not differ strongly 
we expect both in + 2t20 = t and Vn + 2V20 = r to be 
small. 
Expanding to lowest order in t/U we have 



E-{^} 



fk 



u 



(21) 



and if iio is sufficiently small we expect the maximum of 
the lower band to be near (f , f ). At half filling the GS 
energy is 



k 



(22) 



Here a € {10,11,20...} labels the different shells of 
neighbors of a given site and J^ = At'^/U. The 
true expectation value of the Heisenberg antiferromag- 
net in the state |^o) would be obtained by replacing 
—4 (xa ~ j) ~^ Xa — J- For nearest neighbor hopping 
only and xi ~ —0.33 the above result therefore is a fac- 
tor of 2.3 too large. The discrepancy is due to the ne- 
glect of overlap integrals such as dU, see Appendix B. 
The inclusion of such overlap integrals probably is only a 
technical problem but we have to bear in mind that the 
simplification of neglecting such overlap integrals results 
in inaccuracies. 
Next we consider the quasiparticle dispersion relation. 



As already mentioned if tio is small we expect the max- 
imum of the dispersion of the lower band near {^, ^). 
Setting ka = -^ + Ka and expanding to second order in k 
the constant energy contours are elliptical 



E^{k)=E„ 



(k+ - Ko)' 



62 



(23) 



with K± — -j=(kx ± Ky) and - neglecting terms oc t, f as 
well as terms of higher order in t/U - one finds 



1 

n2 



2-^10 (2 ~ 2xio) 
-4tn(^ + 2xii) 



V2 Uio (^ + 2xio) + Jio^-^^il ' ^^'°^4 ~ ^^"^ 



(24) 



Depending on the sign of xio + j the pockets thus are 
shifted towards F or away from F. 

It turns out that the correlation length ( has negligi- 
ble infiuence on the dispersion. More significant is the 
value of the nearest neighbor spin correlation function 
Xio- Figure [1] shows the quasiparticle dispersion and 
the spectral weight of the band for two values of xio- 
Figure [5] shows the pockets obtained for the two different 
values of c/ii 10, which obviously determines the position 
of the hole pocket in the Brillouin zone. The pockets 
are not centered at (f , f ), for the larger values of xio 
the pocket is shifted towards (0, 0) as it is seen in the 
ARPES experiment [1 2] . The overall shape of the dis- 
persion relation is quite similar as the one for a hole in 
an antiferromagnet |42h46| | but it should be noted that 
the present theory does not use any antiferromagnetic 
order. The 'antiferromagnetic shape' of the dispersion 
is due to the fact that even short range antiferromag- 
netic correlations combined with the strong Coulomb re- 
pulsion between electrons are sufficient to suppress the 
nearest neighbor hopping almost completely. More de- 
tailed calculations show that antiferromagnetic spin cor- 
relations instead enhance the incoherent nearest neighbor 
hopping, i.e. nearest neighbor hopping involves emisssion 
or absorption of a spin excitation[47|. It also should be 
noted that the spectral weight of the 'backside' of the 
pocket, i.e. the part facing (7r,7r) is low. This is also 
seen in ARPES |10l - [l3 | ■ although the difference of spec- 
tral weight is much more pronounced there. By analogy 
with hole motion in an antiferromagnet one may assume, 
however, that the coupling to spin excitation and for- 
mation of spin polarons will enhance the difference in 
spectral weight and thus make the theory more similar 
to experiment [48]. 

It should be noted that the antiferromagnetic correla- 
tions determine the location of the pocket in the Brillouin 
zone but not the volume of the Fermi surface, which is 
given by (J18p . As already mentioned, if one sets the spin 
correlation function X — - which actually violates the 
singlet condition ^ - the dispersion relation agrees with 
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FIG. 1: (Color online) Lower panel: dispersion of the lower 
band at 5 = 0.1 for different values of xio- Upper panel: 
dispersion of the spectral weight of the quasiparticle band. 
The correlation length C = 4. 
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FIG. 2: (Color online) Fermi surface at 5 = 0.1 for different 
values of xio- The correlation length C, = A. 



that of the Hubbard I approximation and the hole pocket 
is then centered at (tt, tt). A hole pocket around (tt, tt) has 
indeed been observed in Quantum Monte Carlo (QMC) 
simulations of the Hubbard model]49ll and it is plausible 
that the high temperature used in the (QMC) simula- 
tion renders the spin correlation function small or zero 
and thus shifts the pockets to the corner of the Brillouin 
zone. 

Finally we note that the present theory tends to overes- 
timate the impact of the tn- and t20-terms on the quasi- 
particle dispersion which is why we used relatively small 
values of in and i20 used here. This is probably related 
the fact that no coupling to spin excitations is taken into 
account in the present theory which increases the quasi- 
particle weight and the effect of the in- and i20-terms. 
To conclude this section we return to the issue of the 
hard-core constraint between the holes/double occupan- 
cies. To that end we consider the total densities of 
holes/double occupancies per spin direction: 



nd 



nt 



i 
k 

AT / J ^ ^i*^ l.(T I 



^EK+"^ 



N 



(l-/(i?_(k)))) 



(25) 



This may serve as a criterion for the quality of the ap- 
proximation to relax the hard-core constraint. Namely 
the probability for violation of the constraint at a given 
site is 



p„ = 'indUh 



(26) 



and this is shown in Figure [3] for U /t = 8 as a func- 
tion of the hole concentration 5. This implies that even 
at (5 = 0.2 the constraint is violated at 5% of the sites. 
Enforcement of the constraint e.g. by Gutzwiller projec- 
tion would therefore have a small influence on the results. 
The neglect of the constraint therefore is probably a quite 
reasonable approximation. 



III. SUMMARY AND DISCUSSION 

In summary a theory for the hghtly doped Mott insu- 
lator has been derived. The basic assumption thereby is 
that holes introduced by doping have the same nature as 
the hole-like charge fluctuations at half-filling whose den- 
sity is ex [jyi)- For low hole doping this is probably a rea- 
sonable assumption. The quasiholes then form a Fermi 
gas with a total Fermi surface volume of 5/2. Antifer- 
romagnetic spin correlations render the nearest neighbor 
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FIG. 3: (Color online) Probability p„ for violation of the hard- 
core constraint as a function of the hole concentration S. The 
finite value of p„ at 5 = is due to the charge fluctuations at 
half-filling. 



ture of the phase transition between this correlation- 
dominated low doping phase with a hole-pocket-like 
Fermi surface and the intermediate and electron-density 
phase with a 'large' Fermi surface which has been in- 
ferred e.g. from the dynamical spin and density corre- 
lation function 37| . Experimental data suggest that this 
phase transition occurs at optimal doping or in the over- 
doped range of hole concentrations and thus is related 
to the mechanism of superconductivity. In the frame- 
work of the present formalism this might correspond to a 
replacement of a spin- liquid like 'spin-background' |^o) 
to e.g. a Gutzwiller-projected Fermi sea where the spin 
correlation function has long-ranged Friedel-like oscilla- 
tions. These would introduce long-ranged overlap inte- 
grals between the quasiparticles and thus enhance their 
interaction. 
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hopping essentially ineffective so that the dispersion re- 
lation is similar to the one for hole motion in an antifer- 
romagnet even in the complete absence of static antifer- 
romagnetic order. The Fermi surface thus takes the form 
of four elliptical hole pockets centered near (§,f). In 
addition to these Fermionic excitations the doped insu- 
lator probably has a second type of excitations, namely 
Bosonic spin triplet excitations which are similar in char- 
acter as the magnons at half-filling. We postpone the 
discussion of these excitations and their interaction with 
the charge fluctuations [4 7]. 

A number of simplifying assumptions of different qual- 
ity were made: the spin correlation function of the 'spin 
background' was assumed to have a simple form and was 
taken as a given input parameter. It would be desir- 
able to calculate this e.g. be minimization of the total 
energy but this would require a solution of the system 
of interacting quasiparticles and magnons. On the other 
hand the form of the spin correlation function which was 
assumed seems quite physical and moreover the results 
do not change strongly with the parameters of the spin 
correlation funtion. For example the quasiparticle dis- 
persion is essentially independent of the spin correlation 
legth C- 

In solving the Hamiltonian for the charge fluctuations 
the overlap between pairs of particles has been neglected. 
This is probably the most drastic approximation made 
and was shown to lead to inaccuracies in the results e.g. 
a deviation from the ground state energy at half-filling 
from the known energy of the Heisenberg antiferromag- 
net. The neglected overlaps - being 'four particle over- 
laps' - would create an interaction between the quasipar- 
ticles. The neglect of the hard-core constraint between 
the particles, on the other hand, is probably a very rea- 
sonable approximation because the density of the charge 
fiuctuations is small. 

An interesting question and possibly the key to un- 
derstand high-temperature superconductivity is the na- 



IV. APPENDIX A 

In this Appendix a re-derivation of linear spin wave for 
the spin-i Heisenberg antifcrromagnet is given to show 
the analogy with the present theory for the Hubbard 
model and to some extent justify the neglect of the hard- 
core constraint. For spin wave theory the role of |^o) is 
played by the Neel state, I^^at), and the misalligned spins 
or magnons play the same role as the charge fluctuations 
in the Hubbard model. The misalligned spins are repre- 
sented by Boson operators a^ and bj which are defined on 
the t" and 4,-sublattices, respectively. The a] and b] must 
be chosen as Bosons because the operations of inverting 
spins at different sites commute. These anticommutation 
relations do not hold for operators referring to the same 
site i - rather, for the spin-i system the a'^ and b^ Bosons 
have to obey a hard-core constraint because a spin can 
be flipped only once. Each misalligned spin increases the 
energy by ^ whence we have the representation of the 
longitudinal part: 



Ho 






,ieA 



where z is the number of nearest neighbors and J the 
exchange constant. The transverse part of the Heisenberg 
exchange creates or annihilates pairs of spin fluctuations: 



Hi^ 



J 



E(«^ 



H.c. 



where the sum is over all pairs of nearest neighbors. 
Adding the two terms gives the familiar spin- wave Hamil- 
tonian with a^ and b^ still having to obey a hard-core- 
constraint. This derivation is completely analogous as for 
the charge fluctuations in the Hubbard model. In linear 



spin wave theory the hard-core constraint between the 
Bosons is now simply ignored and the a^ and b^ opera- 
tors are treated as free Boson operators. Despite this, 
Hnear spin wave theory is a highly successful theory and 
the reason is that the density n of Bosons in the ground 
state - obtained self-consistently from the solution of the 
spin wave Hamiltonian itself - is low. Even for the 2D 
Heisenberg antiferromagnet one has n = 0.197 so that 
the probability that two Bosons occupy the same site 
and thus violate the hard core constraint is only n^ « 4%. 
Relaxing the constraint thus will be a very good approxi- 
mation. In the limit of large U and low doping the density 
of charge fluctuations will be small as well (see Figure [3]) 
and we expect that relaxing the hard-core constraint for 
the Fermions will be a reasonable approximation as well. 



V. APPENDIX B 

In this Appendix we show that by properly taking 
into account the overlap integrals the correct expectation 
value of the energy of the Heisenberg antiferromagnet can 



be obtained. We consider half-filling and start with the 
state l^o)- We choose two sites, i and j which are con- 
nected by the hopping term. By acting with the pair 
creation part ex Vij for this bond term we can generate 

the states |1) = dUh] AO) and \2) = d} .h] JO) More pre- 



cisely, the pair creation part generates the state |1) + |2). 
Using the overlap integrals in ([6]) it is straightforward to 
see that this state is an eigenstate of the overlap matrix 
Nij = {i\j) with eigenvalue (1 — Axij) (there is a factor of 
4 due to the prefactor in ([5])) and therefore has the norm 
n — 2(1 — 4:Xij)- The matrix element between j^I'o) and 
the normalized state {l/y/n){\l) + |2)) then is 



1 - 4x, 



so that 2"'* order pertubation theory gives the energy per 
bond Jijixij ~ i)- Here the additional factor of 2 comes 
from the analogous process where the double occupancy 
is created at j and the hole at i. 
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